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ABSTRACT 

Fracture  mechanics  has  traditionally  meant  characterizing  the  linearly  elastic  response 
through  stress  intensity  factors.  However.  It  Is  recognized  that  fracture  depends  strongly 
on  non-linear  material  behavior  in  a  small  region  about  the  crack  tip.  New  experimental 
and  analytical  investigations  must  make  an  effort  to  characterize  this  nonlinear  behavior. 

With  this  observation  in  mind,  we  are  attempting  to  experimentally  investigate  the 
deformation  field  (and  then  calculate  the  strain  and  stress  fields)  at  the  tip  of  a  crack.  The 
experimental  technique,  Moire  Interferometry,  measures  In-plane  displacements  and  Is  sen¬ 
sitive  to  displacements  on  the  order  of  microns.  We  are  now  constructing  an  interferome¬ 
ter  which  will  measure  both  of  the  two  linearly  independent  planar  displacement  com¬ 
ponents  simultaneously;  the  out  of  plane  deformation  can  also  be  determined  on  the  same 
specimen  using  a  Michelson  Interferometer  wrth  similar  sensitivity. 

Much  of  the  data  acquisition  and  reduction  will  be  carried  out  by  computer  through 
digital  image  processing.  The  computer  can  be  used  to  compare  measured  behavior  with 
analytical  or  numerical  predictions  and  help  to  improve  such  models.  We  have  purchased 
the  Image  processing  system  and  are  integrating  the  components.  It  will  be  necessary  to 
write  specialized  software  to  do  the  planned  data  reduction,  but  automating  the  data  pro¬ 
cessing  will  eventually  enable  us  to  examine  more  results  more  thoroughly  than  could  be 
done  by  hand. 

We  have  also  investigated  the  nonlinear  material  behavior  (strain  softening)  numerically 
using  two  different  models.  The  first  is  a  model  of  a  double  cantilever  beam  crack  speci¬ 
men  as  a  beam  on  a  nonlinear  foundation.  The  second  is  a  model  for  the  behavior  of  non¬ 
linear  crazing.  Both  of  these  models  have  been  checked  out  on  the  two-dimensional  prob¬ 
lem  and  were  the  precursor  to  the  problem  of  crack  growth  resulting  front  strain  softening 
material  behavior  which  follows  classically  plastic  behavior. 


1.  INTRODUCTION 


The  fact  that  a  crack  exists  In  a  body  dominates  the  nature  of  the  linearly  elastic  stress 
distribution  near  the  tip  of  that  crack.  This  notion  of  "crack  lip  autonomy"  (with  Its  univer¬ 
sal  inverse  square  root  singularity)  has  traditionally  left  only  the  calculation  of  the 
strength  of  the  singularity  (Stress  Intensity  Factor)  as  the  main  effort  of  fracture  mechanics 
research.  The  stress  intensity  factor  is  a  function  only  of  the  geometry  and  magnitude  of 
the  applied  loading,  and  is  independent  of  material  properties  (for  an  isotropic  linearly 
elastic  solid). 

The  stress  intensity  factor  does  describe  the  conditions  near  the  crack  tip  fairly  well  for 
many  applications,  such  as  fatigue.  However,  it  has  been  recognized  that  the  stress 
response  at  the  tip  of  the  crack  must  be  nonlinear  for  all  but  the  most  brittle  of  materials, 
since  the  stresses  predicted  within  the  linearized  theory  of  elasticity  are  singular,  the 
material  at  the  crack  tip  must  act  in  some  way  to  mitigate  the  effect  of  the  stress  singular¬ 
ity  through  plastic  deformation. 

There  are  several  analytical  models  of  the  nonlinear  (plastic)  behavior  at  the  crack  tip. 
from  the  simple  Dugdale-Barenblatt  cohesive  zone  model  to  the  power  law  hardening 
model  of  Hutchinson,  Rice  and  Rosengren.  Along  with  high  speed  computers  have  come 
numerous  finite  element  and  other  numerical  approximations.  These  models  make 
assumptions  and  predictions  about  the  nature  of  materials  and  their  fracture  response 
which  must  be  tested  under  real  conditions.  In  order  to  test  the  validity  of  the  models, 
actual  displacement,  strain  and  stress  distributions  must  be  determined  experimentally 
and  compared  with  the  predictions.  It  is  also  important  to  determine  the  extent  to  which  a 
basic  assumption,  that  the  material  is  a  continuum,  is  valid;  fracture,  at  its  core,  is  an  inter¬ 
atomic  process  and  hence  a  discrete  process.  Finite  element  models  attempt  to  investigate 
the  crack  tip  region  by  resorting  to  smaller  and  smaller  elements,  each  having  the  proper¬ 
ties  of  the  continuum  material,  but  it  is  obvious  that  there  is  a  scale  beyotid  which  this 
discretization  cannot  be  justified. 

With  these  goals  in  mind,  we  are  attempting  to  measure  with  high  sensitivity  the  in¬ 
plane  as  well  as  the  out-of-plane  displacement  field  in  the  region  of  the  crack  lip.  Once 
the  displacement  field  is  known,  the  strain  field  can  be  calculated  fiom  the  displacement 
gradients  and  the  stress  field  calculated  from  an  assumed  constitutive  law.  These 
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measurements  can  be  used  both  to  validate  numerical  predictions  and  as  a  source  for 
actual  physical  input  in  order  to  improve  the  current  numerical  models.  Our  measure¬ 
ments  will  be  performed  through  Post’s  method  called  Moire  Interferometry,  which  meas¬ 
ures  relative  in-plane  displacements  over  an  entire  field.  The  displacements  are  revealed 
as  fringes,  contours  of  constant  displacement  relative  to  any  other  point  In  the  field. 

Since  each  measurement  will  generate  an  entire  field  of  displacement  information,  it  is 
necessary  to  process  this  data  efficiently.  A  video  digitizing  system  will  thus  be  used  to 
convert  the  fringe  fields  into  a  position-intensity  array  within  a  computer.  The  computer 
will  then  convert  the  fringe  information  to  a  position-displacement  array,  which  can  be 
compared  with  numerical  or  analytical  preditions.  The  computer  can  also  calculate  the  dis¬ 
placement  gradients  numerically  and  create  a  position-strain  array.  With  an  assumed  con¬ 
stitutive  law  the  strains  can  be  converted  to  stresses. 

2.  EXPERIMENTAL  PROCEDURES 

Moire  Interferometry  is  a  highly  sensitive  experimental  technique  which  measures  In- 
plane  relative  displacements  at  the  surface  of  a  specimen.  These  displacements  are 
Inferred  through  the  change  in  the  diffraction  angle  of  a  laser  beam  as  a  diffraction  grating 
deforms.  The  diffraction  grating  is  a  high  frequency  periodic  profile  variation  of  a 
reflective  surface.  The  diffraction  angle  of  the  laser  beam  depends  on  the  wavelength  of 
the  grating  and,  as  the  grating  deforms,  the  grating  wavelength  changes  and  so  does  the 
diffraction  angle.  By  using  two  laser  beams  from  opposite  directions,  very  small  changes 
in  the  diffraction  angles  can  be  detected  by  observing  the  fringe  pattern  resulting  from  the 
Interference  the  two  diffracted  laser  beams.  In  fact,  it  can  be  shown  tliat  each  of  the 
interference  fringes  corresponds  to  a  contour  of  constant  relative  displacement  of  the 
diffraction  grating.  The  displacement  increment  between  adjacent  fringes  is  half  the  initial 
(undeformed)  wavelength  of  the  diffraction  grating.  Tl)e  grating  we  plan  to  use  has  a 
wavelength  of  about  1.7  microns. 

The  required  diffraction  gratings  are  originally  produced  on  high  lesolution  photo¬ 
graphic  plates  such  as  those  used  in  holography.  A  thin  coating  of  aluminum  is  evaporated 
onto  the  surface  of  tlie  plate  to  increase  the  reflectivity.  This  renective  coating  and  the 


grating  profile  can  then  be  transferred  to  any  flat  surface  with  an  adhesive  such  as  epoxy. 
We  have  produced  compact  tension  fracture  specimens  in  this  way,  and  this  method  can  be 
used  on  flat  plate  specimens  of  any  size  and  composition.  Alternatively,  it  is  possible  to 
apply  a  photoresist  to  the  specimen  itself  and  produce  a  grating  directly  on  the  specimen 
in  the  same  way  as  on  a  photographic  plate.  This  should  produce  a  better  quality  grating, 
but  it  limits  the  size  and  composition  of  the  specimen  since  the  photoresist  must  be 
developed  in  a  chemical  bath. 

It  is  necessary  to  hold  two  laser  beams  fixed  in  pace  for  each  component  of  displace¬ 
ment  to  be  measured.  We  are  constructing  an  interferometer  which  can  measure  two 
independent  in-plane  displacements  simultaneously.  As  with  any  interferometric  measure¬ 
ment  system,  vibrations  introduce  noise  to  the  fringe  pattern;  we  hope  to  eliminate  this 
problem  by  mounting  the  interferometer  and  specimen  loading  device  on  a  vibration  iso¬ 
lated  optical  table. 

The  interferometer  itself  uses  three  laser  beams  which  intersect  at  the  surface  of  the 
specimen.  Each  beam  comes  from  one  corner  of  a  square  which  is  parallel  to  the  specimen. 
By  examining  the  interference  between  alternate  pairs  of  intersecting  laser  beams,  three 
different  in-plane  displacement  components  can  be  measured;  0,  45,  and  90  degrees.  (Of 
course  one  of  these  components  can  be  derived  from  the  other  two.)  Each  laser  beam 
must  originate  from  the  same  source  and  all  must  travel  paths  of  similar  length  to  the 
specimen  surface  in  order  to  obtain  high  contrast  interference  patterns. 

Since  each  fringe  represents  a  contour  of  constant  displacement,  as  displacement  gra¬ 
dients  become  high,  so  does  the  fringe  density.  This  could  create  a  problem,  particularly 
for  specimens  exhibiting  highly  plastic  deformations.  However,  with  good  photographic 
equipment  we  expect  to  be  able  to  resolve  100  fringes/millimeter  which  corresponds  to  a 
strain  of  about  1%  for  the  grating  already  described.  This  is  most  likely  adequate  for  most 
metals.  However,  a  coarser  grating  can  allow  for  higher  strain  levels,  and  the  displacement 
sensitivity  can  be  reduced  through  holographic  techniques,  if  that  need  arises  very  close  to 
the  crack  tip. 

At  this  time  the  machinery  for  the  production  of  gratings  on  holographic  plates  is  in 
place,  as  is  the  vacuum  chamber  for  applying  the  reflective  aluminum  coatings.  We  are  able 
to  produce  satisfactory  diffraction  gratings  in  this  manner.  The  process  of  tiansferring  the 
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gratings  to  a  specimen  with  an  adhesive  has  been  tested  and  seems  workable,  the  major 
difficulty  being  air  bubbles  In  the  adhesive  causing  local  faults  In  the  bond.  All  of  the  com¬ 
ponents  of  the  three  beam  Interferometer  are  here  and  the  table  top  loading  device  is  now 
available.  It  will  require  some  time  to  assemble  and  align  the  interferometer  since  precise 
alignment  is  the  most  important  part  of  any  interferometric  measurement  system. 

We  expect  to  have  the  interferometer  set  up  and  some  preliminary  photographs  of 
fringe  patterns  within  the  next  six  months.  We  can  then  decide  If  we  need  to  explore  alter¬ 
native  specimen  geometries  or  different  grating  preparation  methods.  Then  we  can  begin 
to  characterize  the  displacement  fields  as  a  function  of  load,  material  propei  ties  and  crack 
geometry. 


3.  COMPUTER  SYSTEM 

We  have  purchased  a  computer  to  automate  the  reduction  of  the  data  contained  in  the 
interference  fringes.  The  computer  itself  is  a  Digital  Equipment  Corporation  LSI  1 1/73 
central  processing  unit  with  4  megabytes  (MB)  of  random  access  memory  and  an  additional 
180  MB  of  disk  storage  space.  We  also  have  a  70  MB  tape  cartridge  drive  for  permanent 
storage. 


The  fringe  patterns  are  Input  to  the  computer  system  by  digitizing  a  video  image:  the 
analog  video  signal  Is  converted  to  a  digital  representation  by  an  Imaging  Technology 
Incorporated  Analog  Processor  and  sent  along  the  high  speed  video  bus  to  the  Imaging 
Technology  Incorporated  Frame  Buffer  where  it  is  stored  as  an  x-y  position  with  a 
corresponding  brightness  level.  We  have  the  ability  to  convert  a  photograph  to  an  array  of 
1024x960  positions  with  256  brightness  levels.  At  this  resolution,  each  picture  occupies 
about  one  megabyte  of  memory  storage  space.  This  system  has  been  assembled  and  some 
basic  software  to  transfer  images  from  the  frame  buffet  to  the  computer  memory  and  vice 
versa  has  been  written. 


The  next  work  item  concerns  the  development  of  software  to  enhance  and  analyze  the 
fringe  patterns.  The  enhancement  consists  mainly  of  filtering  out  noise  and  thresholding 
the  image  to  reduce  the  fringes  to  thin  black  lines  on  a  white  background  instead  a  varia¬ 
tion  of  gray  levels.  Next,  a  program  to  trace  the  position  of  the  fringes  with  respect  to  the 
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crack  lip  will  be  needed.  This  part  of  the  process  will  require  some  human  Interaction  to 
Impose  a  coordinate  system  and  to  assign  a  displacement  value  to  each  fringe  according  to 
the  boundary  conditions  Imposed.  These  programs  should  also  be  avallale  to  us  within  the 
next  six  months.  The  displacement  field  calculated  can  then  be  compared  to  an  analyti¬ 
cally  generated  displacement  field  at  each  point  in  the  image  array. 

Calculation  of  the  strain  field  in  two  or  three-dimensions  is  straight  forward  in  principle 
once  the  displacement  field  Is  known;  however,  differentiation  can  be  troublesome  as  It 
amplifies  non-smooth  features  In  the  field.  On  the  other  hand,  this  amplification  may  Indi¬ 
cate  regions  where  continuum  as  constitutive  law,  the  stress  field  may  be  calculated.  We 
shall  first  consider  classical  plasticity  theories  such  as  power-lav.'  strain  hardening  models 
to  determine  the  stress  fields.  Imposing  equilibrium  will  also  help  indicate  whether  the 
constitutive  assumptions  are  valid. 

Although  the  initial  emphasis  will  be  on  the  two-dimensional  behavior  of  the  field  quan¬ 
tities.  we  expect  to  proceed  to  the  study  of  three-dimensional  effects  In  a  similar  way. 
These  effects  depend  strongly  on  the  thickness  of  the  specimen  as  well  as  the  similarity  of 
the  crack  geometry  through  the  thickness.  Numerical  calculations  allowing  for  three- 
dimensional  effects  are  also  being  undertaken. 

4.  ANALYTICAL  EFFORT 

4.1.  Background. 

Analytical  studies  in  elastic-jjlastic  fracture  mechanics  have  been  mainly  directed  at 
identifying  certain  fracture  parameters  such  as  tiie  J-integral.  the  crack  tip  opening  dis¬ 
placement  or  angle.  Although  this  approach  has  been  widely  adopted  and  applied  with 
some  success  in  predicting  the  behavior  of  stable  crack  growth  (e.g.  using  the  J-resistance 
curve),  it  fails  to  provide  a  framework  for  studying  in  further  depth  the  effects  of  the 
failure  mechanisms  in  the  crack  lip  zone  on  the  global  behavior  of  crack  growth. 
Specifically,  we  do  not  fully  understand  how  damage-induced  softening  befiavior  of  the 
crack  tip  material  -  which  varies  in  characteristics  with  the  mic  rostruc  tural  f.iilure  mechan¬ 
isms  -  affects  the  stability  of  the  crack  itself.  In  metals,  cracks  advance  with  varying 
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degrees  of  ductility  depending  on  the  loading  and  on  the  environmental  conditions.  If  the 
fracture  is  relatively  brittle,  the  crack  surfaces  separate  easily  and  little  plastic  deforma¬ 
tion  is  incurred  in  the  material  surrounding  the  crack  path.  On  the  other  hand,  when  the 
crack  faces  are  formed  through  void  nucleation,  growth  and  coalescence,  the  effective  sur¬ 
face  energy  turns  out  to  be  substantially  higher  than  the  brittle  surface  energy;  also  the 
extent  of  plastic  deformation  in  the  material  adjacent  to  the  path  traversed  by  the  crack  is 
much  greater.  Thus,  the  material  toughness  and  separation  characteristics  reflecting  the 
failure  mechanisms  at  the  crack  tip  directly  control  the  energy  expenditure  and  the  stress 
and  strain  fields  around  the  crack  front.  Seen  In  this  light,  it  is  clear  that  fracture  analyses 
are  necessarily  Incomplete  without  Incorporating  the  proper  separation  laws  which 
describe  the  failure  modes  at  the  crack  tip;  in  spile  of  the  fact  that  one  parameter  fracture 
criteria  serve  and  have  served  very  useful  purposes  they  seem  to  be  of  somewhat  limited 
virtue  for  understanding  fracture  as  a  fundamental  process.  The  research  we  conducted 
thus  far  has  been  focused  on  these  aspects  of  fracture  discussed  above.  Our  attempts 
have  been  confined  to  planar  geometry  (two-dimensional  problem)  with  extension  to  the 
fully  three-dimensional  case  in  mind. 

Two  pilot  problems,  analyses  of  the  fully  nonlinear  cohesive  zone  model  of  the 
Bareiiblatt-Dugdale  type,  have  been  studied  in  detail.  The  first  problem  is  a  beam  on  non¬ 
linear  foundation  analysis  appropriate  for  a  double  cantilever  beam  (DCB)  specimen  whose 
material  is  elastic  with  the  nonlinear  damage-softened  material  confined  to  a  thin  boun¬ 
dary  layer  on  the  crack  plane  [1].  The  second  problem  concerns  a  nonlinear  craze  embed¬ 
ded  in  an  elastic  medium  subjected  to  symmetrical  loading.  The  nonlinearity  is  again 
confined  to  the  craze  fibril  material  (2|.  These  two  analyses  yield  a  number  of  Interesting 
findings  suggesting  that  the  damage-induced  softening  material  behavior  at  the  crack  tip 
plays  a  central  role  in  determining  the  global  behavior  of  crack  propagation. 

5.  CURRENT  RtSfARCl  I 

We  now  extend  our  study  to  the  case  where  the  material  surrounding  the  crack  is  not 
simply  elastic.  The  first  problem  we  examine  consists  of  a  compact  tension  specitnen 
loaded  by  a  displacement  contiolled  machine.  The  material  is  elastic-plastic  described  by  a 
J2  -How  theory.  The  specimen  is  modeled  by  a  finite  elemetit  mesh  while  the  damage- 


softened  boundary  layer  along  which  the  crack  propagates  is  represented  by  a  row  of  non- 
llnearly  softening  springs.  The  springs'  initial  stiffness  Is  simply  the  material  Young’s 
modulus,  and  the  onset  of  softening  occurs  when  a  critical  nodal  force  Is  reached  and 
failure  takes  place  at  a  critical  separation  distance.  The  "surface  energy"  Is  equal  to  the 
area  under  the  spring's  separation  characteristics;  material  toughness  would  Include  the 
energy  dissipated  in  plastic  deformation.  Note  that  the  definition  of  the  problem  Is  com¬ 
plete  and  no  extra  fracture  criterion  is  needed. 

The  finite  element  program  ABAQUS  is  employed  for  this  purpose.  Preliminary  studies 
indicate  that  the  extent  of  crack  growth  is  critically  dependent  on  the  critical  (threshold) 
nodal  forces  and  on  the  intial  "softening"  slope  of  the  spring. 

More  detailed  studies  are  being  undertaken  and  should  resolve  some  of  the  difficulty 
currently  encountered  such  as  the  proper  scaling  of  the  springs  when  going  from  a  coarser 
mesh  to  a  finer  one.  Early  tests  indicate  that  the  scaling  is  nonlinear  and  may  have  to  be 
obtained  by  trial  and  error.  Finally,  the  nonlinear  spring  technique  can  be  used  in  dynamic 
(not  necessarily  steady  state)  problems  as  well. 
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ABSTRACT 


This  study  investigates  the  effect  of  nonlinear  cohesive  forces  on 
crack  growth  with  the  special  problem  of  craze  mechanics  in  mind.  The 
work  is  presented  in  two  parts.  In  the  first  and  present  one,  we  develop 
a  numerical  method  for  determining  the  equilibrium  shape  of  a  craze  in  an 
infinite  elastic  plane  whose  fibrils  exhibit  very  general  nonlinear 
force-displacement  (P-V)  behavior,  including  strain  softening  characteris¬ 
tics.  The  second  part  of  this  study  [1]^  deals  with  the  numerical  Simula- 

2 

tion  of  craze  and  crack  growth. 


The  problem  formulation  is  based  on  the  superposition  of  the  relevant 
elasticity  Green's  function  and  the  solution  for  the  resulting  nonlinear 
problem  is  effected  by  using  Picard's  successive  approximation  scheme. 
Both  field  equilibrium  and  the  Barenblatt  condition  for  vanishing  stress 
and  strain  singularities  are  satisfied  simultaneously,  rendering  the  craze 
tip  profile  cusp-like.  The  formulation  allows  the  stress  distribution 
profile  and  the  corresponding  P-V  relation  to  be  computed  from  experimen¬ 
tally  measured  craze/crack  displacement  contours  with  certain  advantages 
over  other  methods  proposed  to  date;  it  allows  also  the  computation  of  the 
craze  or  crack/craze  profile  if  the  P-V  relation,  far-fleld  load  and  craze 
or  crack  size  are  specified.  Numerical  investigations  indicate  that  only 
certain  classes  of  the  fibril  P-V  relations  are  consistent  with  measured 
craze  profiles.  In  addition,  it  is  found  that  for  a  given  P-V  relation, 
nontrivial  solutions  exist  only  for  certain  ranges  of  craze  lengths 
depending  on  the  P-V  characteristics  under  consideration. 


1.  The  references  are  listed  at  the  end  of  the  paper. 


2.  The  unabridged  version  of  these  two  papers  is  given  in  reference  2. 
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1.  INTRODUCTION 


When  glassy  polymers  are  subjected  to  extensional  strains  larger  than 
about  0.2\,  tney  begin  to  exhibit  local  'plastic'  deformation.  Two  dis¬ 
tinct  processes  of  plastic  deformations  have  been  observed,  one  being 
mainly  a  shearing  phenomenon  such  as  diffuse  shear  yielding  and  localized 
shear  band  formation,  and  the  other  being  caused  by  cavitation  and  void 
growth  commonly  called  'crazing'.  The  latter  process  is  characterized  by 
a  loss  of  intermolecular  cohesion,  molecular  flow  and  orientation  as  well 
as  a  significant  decrease  in  local  density  [3].  In  this  work,  we  shall 
concern  ourselves  only  with  the  phenomenon  of  crazing. 

The  subject  of  crazing  has  been  extensively  studied  by  both  material 
(polymer)  scientists  and  applied  mechanicians.  Literature  addressing 
various  aspects  of  crazing  can  be  found  in  references  3  through  10.  We 
shall  discuss  here  briefly  the  pertinent  physical  concepts  and  the  motiva¬ 
tions  leading  to  the  present  investigation. 

It  has  been  recognized  for  some  time  that  fracture  in  thermoplastic 
polymers  is  usually  preceded  by  the  formation  of  crazes.  Crazes  are 
formed  in  the  bulk  polymer  along  planes  normal  to  the  direction  of  maximum 
principal  tensile  stress;  they  tend  to  nucleate  at  micro-defects  or  inclu¬ 
sions  and  grow  by  the  formation  of  fibrils  from  the  bulk  polymer  [6,7]  and 
thicken  through  a  combination  of  two  distinct  mechanisms,  namely  through 
'surface  drawing'  and  'creep' .  In  surface  drawing,  new  polymer  is  drawn 
from  the  bulk  into  the  craze  in  the  form  of  fibrils.  During  creep  exist¬ 
ing  fibrils  are  simply  elongated,  mostly  locally  around  the  weak  portions. 

Note  that  although  'creep'  implicitly  signifies  a  time  dependent  process, 
we  will  exclude  time  consideration  from  this  work.  For  our  purpose  it  is 
sufficient  to  recognize  that  the  interplay  of  these  two  mechanisms  gives 
rise  to  an  effective  force-displacement  (P-V)  behavior  which  we  consider 

3.  For  an  analysis  accounting  for  these  two  basic  mechanisms,  see  the  recent 
paper  by  Kramer  and  Hart  [11]. 
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to  be  Initially  prescribed.  Here  'force'  represents  the  cohesive  stress 
exerted  by  the  fibril  on  the  craze/bulk  polymer  boundary.  'Displacement' 
refers  to  the  net  displacement  of  the  craze  boundary  (see  Section  2.1  for 
details) . 

Electron  microscopy  shows  that  the  typical  craze  thickness  is  on  the 

order  of  1  and  the  length  to  thickness  (aspect)  ratio  is  approximately 
2  3 

10  to  10  .  Crazes  generally  consist  of  a  network  of  fibrils  between  10 
to  40  nm  in  diameter,  Interspersed  with  voids  of  similar  dimensions.  The 
transition  boundary  between  a  craze  and  the  bulk  polymer  where  surface 
drawing  takes  place  is  very  thin  i.e.  less  than  2  nm.  Craze  fibrils  have 
been  shown  to  exhibit  considerable  strength  [8,9].  Therefore,  crazing 
could,  under  certain  conditions,  substantially  enhance  the  toughness  of 
bulk  polymers.  However  the  presence  of  crazes  does  not  necessarily 
enhance  the  strength  of  the  overall  structure.  In  many  instances,  uncon¬ 
trolled  craze  growth  tends  to  Induce  fracture,  causing  structures  to  fall 
'prematurely'.  For  more  in-depth  discussions  on  microstructures  and  pro¬ 
perties  of  crazes,  see  references  7  and  10.  Figure  1  shows  a  typical  craze 
and  its  tip  region.^ 

We  next  review  briefly  the  analytical  approaches  that  have  been 
attempted  in  modelling  the  crazing  phenomenon  to  date.  Barenblatt  con¬ 
sidered  the  problem  of  equilibrium  cracks  in  brittle  fracture  [12];  even 
though  the  analysis  was  carried  out  in  the  context  of  cracks,  the  basic 
propositions  and  concepts  apply  to  crazes  equally  well.  His  proposition 
eliminates  the  infinite  stresses  and  strains  at  the  crack  tip  by  incor¬ 
porating  cohesive  stresses  distributed  over  a  cohesive  zone  of  finite 
length  such  that  the  crack  is  in  "equilibrium"  and  the  stress  Intensity 
factor  vanishes.  Assuming  a  constant  cohesive  (yield)  stress  in  this 
cohesive  zone,  Dugdale  derived  a  simple  relationship  between  its  size  and 
the  applied  stress  [13]  in  plane  stress;  this  study  was  then  considered  by 

4.  The  micrographs  were  taken  by  L.  Berger  and  made  available  to  us  by  Pro¬ 
fessor  E.J.  Kramer. 


Rice  In  more  detail  using  the  rigid-plastic  strip  model  [14].  Goodler  and 
coworkers  calculated  the  rate  of  plastic  work  dissipation  for  the  same 
model  with  an  added  assumption  that  a  crack  propagates  with  a  self-similar 
geometry  [15]  and  later  introduced  a  crack  propagation  model  using  non¬ 
linear  atomic  separation  laws  indicating  how  the  cohesive  stress  distribu¬ 
tion  changed  as  the  crack  extended  [16].  Atkinson  proposed  an  iterative 
scheme  for  solving  axially  symmetric  problems  relating  to  cracks  that  open 
under  a  displacement-dependent  internal  stress  without,  however,  satisfy¬ 
ing  the  Barenblatt  condition  [17].  Andersson  and  Bergkvist  solved  a  crack 
problem  using  a  strip  model  [18]:  Their  strip  material  behavior  was 
plecewlse-llnear  with  a  "linear  softening"  characteristic;  the  strip 
extended  to  infinity  and  the  displacements  decreased  gradually  to  zero,  so 
that  the  stresses  were  bounded  everywhere,  similar  to  those  in  reference 
16.  E.  Smith  derived  a  class  of  acceptable  force-displacement  relation¬ 
ships  by  starting  with  displacement  profiles  chosen  so  as  to  yield  simple 
solutions  to  the  integral  equation  which  relate  the  displacement  and 
stress  within  the  cohesive  zone  [19]  under  satisfaction  of  Barenblatt 's 
boundedness  condition. 

In  addition  to  these  stress  analyses,  there  are  a  number  of  contribu¬ 
tions  concerned  with  the  mechanics  of  craze  growth.  In  1965,  Knight  [20] 
employed  the  Fourier  transform  method  [21]  to  calculate  the  stress  distri¬ 
bution  along  a  craze,  but  due  to  the  erroneous  craze  profile  used,  the 
resulting  stress  distribution  was  incorrect.  Knauss  presented  a 
boundary-layer  model  (formulated  for  a  semi-infinite  crack)  which  is 
extendable  to  account  for  rate-dependent  cohesive  stress-displacement 
behavior  [22].  The  critical  strain  (displacement)  criterion  was  employed 
and  it  was  found  that  for  materials  with  softening  force-displacement 
characteristics,  the  crack  tip  cohesive  stress  decreased  rapidly  to  zero 
while  the  size  of  the  cohesive  zone  rose  sharply  as  the  critical  strain 
was  approached  at  the  crack  tip. 


Verheulpen-Heymans  and  Bauwens  presented  a  modified  Dugdale  model 
with  a  constant  cohesive  stress  acting  over  a  short  distance  behind  the 
craze  tip  and  a  lower  constant  cohesive  stress  acting  over  the  rest  of  the 
craze  length  [23].  Lauterwasser  and  Kramer  performed  experiments  in  which 
accurate  measurements  of  the  density  of  the  craze  materials  together  with 
the  craze  profile  permitted  the  determination  of  the  stress  distribution 
profile  [7]  through  use  of  the  Fourier  transform  technique  [21].  They 
also  demonstrated  conclusively  that  surface  drawing  was  the  dominant  pro¬ 
cess  in  craze  thickening.  As  we  shall  see  later,  this  finding  is  of  gen¬ 
eral  significance  in  determining  the  appropriate  P-V  relations  for  craze 
fibrils.  Subsequently,  Wang  and  Kramer  employed  a  distributed  dislocation 
model  to  calculate  the  far-fleld  applied  stress  as  well  as  the  surface 
stress  profile  along  the  craze  from  the  experimental  displacement  profile 
in  the  craze  zone  [24].  The  method  gave  results  in  good  agreement  with 
the  Fourier  transform  method.  The  stress  profiles  obtained  in  this  way 
for  polystyrene  and  polycarbonate  exhibited  stress  concentrations  at  the 
crack  tip  as  well  as  at  the  craze  tip  and  approximately  constant  stress  in 
between.  Donald  and  Kramer  later  examined  four  additional  polymers  and 
obtained  similar  results  [25]. 

Attempts  to  determine  stress  distribution  profiles  for  given  dis¬ 
placement  profiles  by  using  other  methods  such  as  finite  element  and  boun¬ 
dary  element  analyses  were  made  by  Bevan  [26]  and  Sun  and  Hsiao  [27], 
where  the  latter  authors  also  allowed  yielding  in  the  bulk  polymer  in 
their  finite  element  study.  Although  the  results  of  both  studies  agreed 
generally  with  those  obtained  by  Kramer  et  al  [7]  through  the  Fourier 
transform  method,  the  accuracy  and  efficiency  seems  wanting. 

So  far  we  have  reviewed  studies  in  which  the  stress  distribution 
along  the  craze  boundary  is  calculated  from  the  experimental  displacement 
profile  by  using  various  analytical  techniques.  The  main  purpose  of  the 
present  study  is,  however,  to  address  the  more  difficult  problem  where  one 
seeks  solutions  for  both  the  stress  and  displacement  profiles  together 
with  the  far-fleld  applied  stress,  given  a  craze  length  and  a  nonlinear 


fibril  P-V  relation,  such  that  Barenblatt's  condition  Is  also  satisfied. 
This  problem  has  thus  far  not  been  solved  correctly.  Recently,  Walton  and 
Weitsman  presented  a  solution  scheme  for  the  special  case  In  which  the 
fibril  P-V  relation  Is  linear  [28]  and  later  Weitsman  [29]  extended  that 
analysis  to  Include  the  case  where  the  P-V  relation  Is  nonlinear.  The 
analyses  in  references  28  and  29  are,  however,  deficient  In  several 
aspects.  For  instance,  their  solution  did  not  satisfy  the  Barenblatt  con¬ 
dition  but  had  to  be  'corrected'  later  by  imposing  additional  cohesive 
stresses  in  a  small  'tip  zone'  to  make  the  stress  Intensity  vanish.  This 
correction  Is  artificial  and  results  In  stress  distributions  that  are 
inconsistent  with  those  deduced  from  experiments  [24,25].  This  shortcom¬ 
ing  Is  due  partly  to  an  Inaproprlate  choice  for  the  fibril  P-V  behavior; 
Interested  readers  may  wish  to  consult  reference  2  for  a  more  complete 
discussion. 

In  view  of  these  past  contributions,  it  becomes  clear  that  an  analyt¬ 
ical  and  computationally  efficient  method  for  determining  the  displacement 
and  stress  profiles  of  a  craze  with  prescribed  nonlinear  fibril  behavior 
under  satisfaction  of  both  equilibrium  and  the  smooth-closure  conditions 
is  still  lacking.  The  primary  objective  of  this  work  is  to  present  such  an 
analysis  and  to  use  it  in  the  study  of  various  aspects  of  the  mechanics  of 
craze  and  crack  growth  in  thermoplastic  materials. 

In  the  following  sections,  we  discuss  the  model  for  representing  a 
craze.  The  idealizations  made  to  model  the  problem  are  explained  in  Sec¬ 
tion  2.1.  The  mathematical  formulation  including  nondimenslonalizatlon 
and  discretization  for  final  implementation  then  follow  in  Sections  2.2 
and  2.3.  In  Section  3,  two  alternate  methods  for  computing  stress  distri¬ 
bution  profiles  from  craze  displacement  contours  are  presented  for  the 
case  of  a  full  craze  and  for  a  craze  with  a  central  crack.  Because  of 
their  simplicity,  these  methods  are  both  computationally  more  efficient 
and  inherently  more  accurate  (for  the  same  degree  of  discretization)  than 
those  used  in  the  past  by  other  investigators.  The  convergence  charac¬ 
teristics  of  the  solution  technique  are  briefly  explored  and  the  general 


applicability  and  performance  of  the  numerical  schemes  are  assessed  In 
Section  4.1.  Some  characteristics  of  the  fibril  P-V  relations  that  are 
inconsistent  with  realistic  craze  profiles  are  then  examined  and  the 
existence  of  solutions  is  addressed  briefly  in  Section  4.2.  The  question 
of  the  uniqueness  of  solutions  is  raised  but  not  addressed  rigorously. 

2.  PROBLEM  FORMULATIONS  AND  SOLUTIONS 

In  this  section,  we  reduce  the  physical  problem  to  a  mathematically 
manageable  one  through  suitable  approximations.  The  resulting  equations 
are  then  nondimenslonalized  and  solved  through  a  numerical  algorithm  based 
on  Picard's  iteration. 

2.1.  Geometric  and  Mathematical  Idealization. 

We  consider  first  the  geometry  of  the  problem  as  depicted  in  Figure 
2a.  A  craze  of  length  C  with  an  (internal)  central  crack  of  length  A  is 
embedded  in  an  infinite  plate  subjected  to  a  remote  and  uniform  tensile 
stress  applied  normal  to  the  major  axis  of  the  craze.  X.Y  are  the 
coordinates  and  W(X)  denotes  the  Y-displacement  of  the  craze  contour.  We 
let  P(X)  represent  the  cohesive  force  (stress)  distribution  exerted  by  the 
fibrils. 

As  shown  by  Lauterwasser  and  Kramer,  a  craze  is  formed  from  a  primor¬ 
dial  craze,  characterized  by  a  contour  W^(X);  this  contour  contains  the 
bulk  polymer  material  which  Is  subsequently  transformed  into  craze  fibrils 

through  the  surface  drawing  process  [7].  Thus  the  dashed  contour  W  (X)  in 

o' 

Figure  2b  is  displaced  in  the  drawing  process  to  the  solid  contour  W(X) 
due  to  the  application  of  Since  we  know  from  experiments  that  Wq(X) 
is  much  smaller  than  C  and  that  the  slope,  W'^(X),  is  also  very  small 
everywhere,  we  can  recast  the  geometry  of  the  problem  as  shown  in  Figure 
2c  where  V(X)  is  the  net  displacement  of  the  craze  profile  caused  by  the 
applied  stress  as  measured  from  W^(X)  so  that  V(X)  =  W(X)  -  W^(X).  Cast 


9 


In  this  light,  the  force-displacement  (P-V)  relation  may  be  represented 
schematically  as  shown  In  Figure  2d.  Keeping  in  mind  that  V  <  0 
corresponds  to  the  primordial  craze  under  small  strain  prior  to  fibrilla¬ 
tion,  the  dashed  portion  of  the  P-V  curve  Is  assumed  linear  with  a  slope 
equal  to  Young's  modulus,  E,  for  the  bulk  polymer.  At  V  =  0,  the  thres¬ 
hold  (or  fibrillation)  stress  P^  Is  reached  and  from  then  on  the  P-V  rela¬ 
tion  Is  nonlinear.  Three  P-V  relations  are  Illustrated,  curve  #1  is 
representative  of  a  craze  that  does  not  contain  a  crack.  We  shall  later 
refer  to  such  a  craze  as  a  'full  craze'.  Curve  #2  allows  for  continued 
hardening  before  softening  to  failure  (fibril  degradation  and  rupture)  at 
V^.  Curve  #3  Is  typical  of  a  craze  with  a  central  crack;  this  rehardening 
behavior  may  be  understood  in  the  context  of  the  midrib  behavior  [7].^ 

2.2.  The  Mathematical  Formulation. 

By  superposing  the  Green's  functions  obtained  through  Westergaard 
stress  function  method  [30],  the  solutions  to  the  problem  can  be  written 
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where  7f  »  (j:^)  for  plane  stress  and  (3-4»')  for  plane-strain;  v  = 

Poisson's  ratio;  n  =  shear  modulus.  is  the  mode-I  stress  intensity 

5.  The  midrib  Is  formed  initially  at  the  craze  tip  due  to  fibril  drawing  at 
high  stress.  As  the  craze  grows,  the  midrib  is  unloaded  as  it  becomes 
part  of  the  mature  (developed)  craze.  The  introduction  of  a  crack  In  the 
craze  causes  the  stress  at  the  crack  tip  to  rise,  so  that  the  midrib  Is 
reloaded  and  therefore  must  sustain  higher  stress. 


6.  Huang  independently  employed  an  Identical  formulation  in  studying  the  size 
of  the  crack  tip  cohesive  zone  based  on  an  Interatomic  potential  [31]. 


factor,  and  V(X)  Is  the  net  displacement  due  to  the  far-fleld  loading  and 
the  fibril  cohesive  stress. 


For  presentation  of  the  results.  It  Is  helpful  to  non-dimenslonallze 

pertinent  parameters.  We  let  a  craze  length  be  the  reference  length. 

The  stresses  are  scaled  by  (-r^) ;  in  this  fashion,  the  Poisson's  ratio  Is 

1+J7 

absorbed  for  both  cases  of  plane  stress  and  plane  strain.  In  short. 


for  plane  stress 
for  plane  strain 


We  shall  use  lower  case  letters  for  dimensionless  quantities  and 
define 


Thus,  in  dimensionless  forms,  (2.1)  and  (2.2)  are 


c 

v(x)  =■  ;  [<t^-p(t)]  m(c,x,T)  dr 
o 


-  r 


i 

where  the  kernel  in  equations  (2.3 


c 

I 


o 

) 


<^«,-P(t) 

— - dr 

is  given  by 


1 

I 

I 

J 


(2.3) 


(2.4) 


m(c,x,T) 


,  I  2  2  .  2  2 

\|c  -X  -t-\  c  -T 


V  2  2.12  2 

\  c  -X  -\  c  -T 


(2.5) 


Note  that  here  k.  is  dimensionless.  Integrating  equation  (2.3)  yields 
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v(x)  =  <T^  \|c^-x^  -  /  m(c,x,T)  p(t)  dT 

o 


(2.6) 


For  the  stress  at  the  craze  tip  to  be  finite,  or  equivalently,  for 
the  craze  tip  to  close  smoothly,  i.e.  v'(c)  =  0,  it  is  necessary  and  suf¬ 
ficient  that  kj  vanish.  Then  (2.4)  gives 


^  ,2^  -J^dT 

oo  7t  - 

°  \ 

Substituting  a  from  (2.7)  into  (2.6).  one  obtains 

OO 

c 

v(x)  =  /  h(c,X.T)  P(t)  dT 
o 


2  2 
c  -T 


(2.7) 


(2.8) 


where 


"(^■*••'1  ■;  1^  -  (2.91 

The  formulation  is  completed  by  specifying  the  fibril  force-displacement 
behavior 


p(x)  -  p[v(x)] 


(2.10)’^ 


Equations  (2.8)  and  (2.10)  must  be  solved  simultaneously.  This  is  accom¬ 
plished  by  using  Picard's  method  of  successive  approximation  outlined 
below® 


7.  The  more  general  case,  p(x)  =  p[v(x),x]  can  be  handled  within  the  frame 
work  of  our  formulation  using  the  same  solution  scheme.  A  p-v  relation 
with  explicit  dependence  on  x  is  a  more  realistic  candidate;  however  due 
to  lack  of  any  such  data,  we  confine  our  attention  to  the  form  (2.10). 


8.  As  it  turns  out,  this  simple  scheme  is  much  more  efficient  than  the  Newton 
iterative  scheme  which  was  also  Implemented  and  yielded  the  same  results. 
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We  start  with  an  initial  guess  of  the  form, 

v°(x)  =  o°(c  -  x)"'  (2.11) 

oo 

where  0°  and  r  are  to  be  appropriately  chosen;  other  forms  of  v°(x)  are 
also  acceptable. 

The  approximation  is  given  by 

v^(x)  =  /  h(c,x,T)p[v'^'^(T)]dT  ;  k  =  1.2 -  (2.12) 

We  consider  that  a  solution  has  been  attained  when  v  (x)  differs  from 
k-1 

V  (x)  at  all  (discretized)  points  by  0.1%  or  less.  The  convergence  of 
the  scheme  is  found  to  be  relatively  insensitive  to  the  choice  of  the  ini¬ 
tial  guess  of  and  r  (see  Section  3.1). 

It  is  worth  noting  that  in  the  case  where  a  crack  is  present  as 
determined  by  the  critical  crack  tip  opening  displacement  (CTOD)  criterion 
(see  [1]  for  further  discussion),  the  crack  length  is  not  fixed  a  priori 
but  is  obtained  as  part  of  the  solution.  We  next  discretize  the  equations 
for  numerical  implementation. 

2.3.  Discretization  and  the  Resulting  Algorithm. 

Because  the  kernels  ra(c,x,T)  and  h(c,x.T)  defined  by  (2.5)  and  (2.9) 

can  be  evaluated  in  closed  forms,  the  Integration  of  these  'singular' 

Integrals  can  be  calculated  very  accurately  without  requiring  excessive 

g 

subdivisions  of  the  domain.  We  discretize  the  domain  of  integration  into 
n  equal  Intervals  of  length  Ax  ,  where 

9.  For  contributions  about  the  singular  points,  the  Cauchy  Principal  values 
of  the  integrals  are  used. 
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Ax  =  AT  =  c/n 

and  are  evaluated  at  the  midpoint  of  each  interval, 

in  mind,  an  integral  I(x)  can  be  approximated  as^® 

c 

I(x^)  S  /  f(c.X^,T)  p(T)dT  =  fj(C.X^)Pj  =  f^jPj 

where 

Tj+At/2 

fi.  S  f  (c,X.)  S  /  f(C,X.,T)  dT 

Tj-At/2 

Integrals  of  the  form  (2.14)  in  this  study  are  easily  evaluated 
forms;  specifically,  from  (2.5)  and  (2.6),  one  obtains 


“l  ■  -  "ijPj 


where 


m. j  H  -  I  T  log 


C^-T^ 


M 


12.  2 
-x^-\  c  -1 


2^|c^-XjCOS~^(^) 


X^log 


Xj\ 


2  2  12  2 
C  -T  -T^jc  -X^ 


Similarly,  (2.7)  may  be  written  as 


T=T  .+At/2 

1  J 


X. \ |c^-T^+T^ jc^-X^  |T=Tj-AT/2 


With  this 

(2.13) 

(2.14) 

in  closed 

(2.15) 


(2.16) 


10.  Repeated  indices  indicate  summation  from  1  to  n. 
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<T^  =  ;  l(C.T)  p(T)dT  =  I  p 
O 


(2.17) 


2  1 

where  1(c,t)  =  -  - - -  and 

71 


\|c2-r2 


Or  -At  T=T.+AT/2 
1  =  -(  gin  ^(-)1  ^ 

j  nL  ®  c'Jt=t  -At/2 


(2.18) 


In  view  of  (2.8)  and  (2.9),  we  obtain  in  compact  form. 


V.  -  V(x.)  H  h^.Pj 


(2.19) 


where 


h  =  ^c^-x^  I 

ij  “  \r  1  j 


■“ij 


(2.20) 
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It  turns  out  that  from  equations  (2.16),  (2.18),  and  (2.20),  T.  is 

independent  of  c  while  m^^  and  h^j  are  homogeneous  in  c.  These  properties 


prove  to  be  very  useful  when  we  employ  the  numerical  algorithm  to  study 
craze  and  crack  growth  in  the  second  part  of  this  work  [1].  In  concise 
form,  one  has 


Ij  =■  lj(c,n)  =  lj{l,n) 


(2.21) 


"iJ  "  ^  *  "ij(l-n) 


(2.22) 


^ij  "  hjj(c,n)  =  c  •  h^jd.n) 


(2.23) 


Thus,  for  a  given  degree  of  discretization  n  ,  d,  and  h..  only  need 

to  be  computed  once  for  c=l  and  stored  for  later  use.  In  the  simulation 


11.  Note  that  the  matrix  h  is  singular  l.e.  given  a  displacement  profile 
v(x),  the  cohesive  stress  profile  can  only  be  determined  up  to  an  additive 
constant.  This,  however,  does  not  affect  the  validity  of  our  algorithm. 


-•/-"/A 


15 


of  craze  and  crack  growth  [1],  c  Increases  continually,  thus  the  proper¬ 
ties  exhibited  in  (2.21-2.23)  reduce  the  computation  time  substantially. 
Note  that  1^,  and  h^j  are  independent  of  the  p-v  relation  considered. 

The  algorithm  can  now  be  summarized.  Given  a  craze  length  c  and  a 
p-v  relation,  we  seek  solutions  in  the  form  of  p{x),  v(x)  and  <J^  according 
to  the  following  procedure: 


Prescribed  p-v  relation  : 

p(x. ; 

1  =  Pi  = 

p(v(x^) 

)  =  p(Vi) 

(2, 

,24a) 

Initial  guess  ; 

v°(x 

\  0 
=  ''i 

i 

-  <I  (c-x 

OO 

(2 

.24b) 

k'‘^-iteration  : 

v'^  = 
1 

v''(X.) 

o 

II 

\  /  k-l, 

n)p(Vj  ) 

(2, 

.24c) 

P(V^) 

(2, 

.24d) 

/  . 
OO 

T  ^ 

1  .p . 

(2, 

,24e) 

Convergence  is  considered  achieved  when  the  maximal  change  for  any  v^  is 
less  than  0.1%  between  two  successive  iterations. 


3.  SCHEMES  FOR  DETERMINING  STRESS  PROFILES  FROM  DISPLACEMENT  PROFILES 

As  mentioned  in  the  Introduction,  for  the  case  of  a  full  craze,  the 
stress  distribution  profile  p(x)  may  be  obtained  from  the  measured  dis¬ 
placement  profile  v(x)  and  using  the  Fourier  transform  method  [21]. 
However,  for  a  craze  with  a  central  crack,  only  the  displacement  profile 
in  the  craze  zone  is  needed  to  determine  both  p(x)  and  using  a  distri¬ 
buted  dislocation  scheme  [24].  Here  we  present  two  alternative  methods 
for  determining  the  stress  profiles  based  on  the  formulation  in  the  previ¬ 
ous  section;  one  method  is  suitable  for  a  full  craze,  the  other  for  a 
craze  containing  a  central  crack.  Both  schemes  have  been  fully  tested  and 
found  to  be  numerically  more  efficient  than  the  methods  previously 


proposed  [7,24,26] 


First  we  discuss  the  method  suitable  for  the  problem  of  a  full  craze. 
Given  v(x^)  and  a^,  equation  (2.15)  can  be  rewritten  as 


—  I  2  2 

m.  .p .  =  (<J  \  |c  -X.  -  V. ) 

'  “"Xj  1  i' 


(3.1) 


Equation  (3.1)  la  a  system  of  n  linear  equations  from  which  p^.  can  be 
determined  either  by  Gaussian  elimination  or  by  inverting  m^j  .  in  the 
latter  case,  (3.1)  and  (2.22)  gives 


Pj  =  m'^(c,n){o_^^|c^-x.  -  v^) 


(3.2) 


where 


mjj(c,n)  =  m"j(l,n)/c 


(3.3) 


Note  that  m~j(l,n)  only  needs  to  be  computed  once  and  stored  for  later 


For  demonstration  purposes,  this  scheme  is  applied  to  v(x)  and  as 
presented  in  reference  7  (see  Figure  3).  The  resulting  cohesive  force 
distribution  p(x)  shown  in  Figure  4  and  that  derived  in  [7]  through  the 
Fourier  Transform  method  are  virtually  identical.  The  p-v  relation 
obtained  by  combining  Figures  3  and  4  is  plotted  in  Figure  5. 


This  method  is,  of  course,  also  applicable  to  a  craze  with  an  inter¬ 
nal  crack.  However,  since  the  measured  displacement  profiles  in  the  crack 
region  are  usually  not  reliable  on  a  point-by-point  basis  (see  [24]),  an 
alternative  scheme  using  only  the  displacement  profile  in  the  craze  zone 
as  boundary  data  is  preferable.  Such  a  method  has  been  presented  by  Kra¬ 
mer  and  Wang  [24];  however,  the  method  requires  differentiation  of  the 
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experimentally  measured  displacement  profile  which  incurs  additional 
error.  In  view  of  this  shortcoming,  we  present  here  a  scheme  which 
requires  as  inputs  only  the  displacement  in  the  craze  zone  and  the 
crack  tip  location.  The  outputs  are  the  stress  profile  p(x) .  the  applied 
stress  <T^  and  the  displacement  in  the  crack  region. 

Let  i^  be  such  that  the  fibril  cohesive  force  at  the  i^^  station  is 


Pi 

=  0 

for 

i  <  i 

Pi 

0 

for 

1  >  i 

Given  the  displacement  in  the  craze  zone,  v^  for  i  t  i^,  we  wish  to  find 
and  p^  for  1^1^  and  v^  for  i  <  1^  .  Consider  equation  (2.19), 

Vi  =  iijjPj  ,  i,j  =  1,2 . n  (2.19) 

In  light  of  (3.4),  (2.19)  reduces  to 

Vj  =  0  :  1  =  1,2,...  ,1^.-1  (3.5a) 

^1  =  ;  l.k  =  i^,i^.  +  l . n  (3.5b) 

Equation  (3.5b)  can  be  solved  for  p^^,  °  . "  Gaussian  elimina¬ 

tion.  Then  by  (2.17)  and  (2.19) 

<^oo  =  ijPj  =  ik^k  ’  ^  "  H . "  (3-®) 

Vi  =  1  1 . ^  "  ^t . " 

This  completes  the  presentation  of  the  two  alternative  schemes  for  comput¬ 
ing  p(x)  from  v(x) . 


4.  THE  CHARACTERISTICS  OF  THE  ALGORITHM  AND  THE  p-v  RELATION 


We  examine  next  two  aspects  of  the  craze  problem  from  a  general 
viewpoint.  First  we  report  on  the  convergence  characteristics  of  the 
algorithm  and  attempt  to  give  some  broad  criteria  on  its  applicability  and 
effectiveness.  Then  two  models  of  p-v  relations  are  studied  with  the  aim 
of  determining  how  the  existence  of  a  solution  depends  on  the  nonlinearity 
of  the  p-v  relation.  The  findings  Indicate  that  there  are  restrictions  on 
the  characters  of  the  p-v  relations  if  realistic  craze  profiles  are  to 
result,  i.e.  profiles  with  nonnegative  (non-interpenetrating)  displace¬ 
ments  at  all  points. 

4.1.  The  Convergence  of  the  Algorithm. 

Having  discussed  the  algorithm  and  the  convergence  criterion  in  Sec¬ 
tions  2.2  and  2.3,  we  consider  the  characteristics  of  the  algorithm  as 
applied  to  specific  problems.  In  Figures  6,  we  show  how  the  solution  for 
v{x)  for  a  (half-cosine)  p-v  relation  converges  for  an  initial  guess  of 
the  form 


V°(X)  =  <T°(C  -  X)" 

where  in  this  example,  the  initial  guesses  are  =  0.05  and  r  =1.5 

OO 

Nine  Iterations  are  required  for  convergence  in  this  case.  (Each  itera¬ 
tion  takes  about  1.2  CPU  seconds  on  a  VAX  11-780.)  The  solution  for  <r  is 

OO 

0.0176.  Again,  a  displacement  tolerance  of  O.lSi  is  used  as  the  convergence 
criterion.  For  n=200,  we  find  that  checking  for  convergence  at  only  20  or 
100  equi-spaced  points  give  identical  results  to  checking  at  all  200 
points.  This  is  so  because  the  convergence  is  quite  uniform  for  most  p-v 
relations.  Throughout  this  work,  we  conservatively  use  100  checkpoints 
for  n=200.^^ 

12.  Specifically,  the  p-v  relation  in  Figure  7(b)  with  p  =  p  =  0.0393, 

V  =0.0  and  v  =  0.008  is  used.  ° 

m  c 
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With  regard  to  the  general  applicability  of  the  algorithm,  we  find 
that  the  initial  guess  does  not  need  to  be  very  close  to  the  correct  solu¬ 
tion  for  the  scheme  to  work  as  illustrated  in  Figure  6.  It  is  found  that 
the  exponent  'r'  in  equation  (2.24b)  for  the  initial  guess  v°(x)  may  range 
from  0.5  to  2.0  but  r=1.5  seems  to  work  best  in  most  instances.  When 
there  is  no  solution  such  as  those  situations  to  be  discussed  later  in 
Section  4.2,  the  algorithm  quickly  diverges  for  a  wide  range  of  initial 
guesses  (see  Section  3.2  of  reference  2  for  details).  For  p-v  relations 
which  soften  initially  and  then  reharden  such  as  the  one  studied  in  Part  2 
of  this  work,  the  convergence  could  be  very  slow  and  30-70  Iterations  may 
be  required  depending  on  the  initial  guess. 

In  cases  where  the  algorithm  did  not  converge,  the  Newton  iterative 
scheme  was  tried  but  did  not  yield  convergence  either.  Per  iteration,  the 
Newton  scheme  takes  about  250  times  the  CPU  time  needed  for  Picard's 
method.  For  all  cases  studied  Picard's  iterative  scheme  yielded  the  same 
results  as  the  Newton  method. 

4.2.  The  Characteristics  of  Admissible  p-v  Relations. 

We  examine  next  how  the  various  characteristics  of  a  given  nonlinear 
p-v  relation  affect  the  existence  of  a  solution.  This  issue  is  of  consid¬ 
erable  Importance  as  little  is  known  regarding  the  fibril  p-v  behavior 
that  is  consistent  with  observed  craze  displacement  profiles.  When 
'improper'  p-v  behaviors  are  used,  the  results  obtained  are  necessarily 
incorrect  as  discussed  in  Section  1  concerning  recent  works  in  references 
28  and  29.  (E.  Smith  [19]  and  Knauss  [22]  presented  some  admissible  p-v 
relations  for  the  crack  tip  cohesive  zone.) 


13.  We  tried  n=50,  100,  200,  300,  400,  500  and  found  n=200  to  be  the  best 

compromise  between  the  Integration  accuracy  (within  0.0001%)  and  the  com¬ 
puting  time. 


V 
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In  the  present  study,  we  find  that  there  are  certain  restrictions  on 
the  p-v  relations  for  the  algorithm  to  converge.  In  reference  2,  heuristic 
arguments,  similar  to  the  one  presented  in  connection  with  the  finding  (1) 
below,  were  used  to  demonstrate  that  the  proposed  algorithm  is  capable  of 
seeking  out  a  solution  for  a  given  p-v  relation  and  craze  length,  if  such 
a  solution  exists.  Otherwise,  our  numerical  studies  show  it  to  be  most 
likely  that  algorithm  divergence  indicates  the  absence  of  a  solution.  The 
question  of  uniqueness  is  then  addressed,  but  only  a  brief  summary  of  the 
findings  is  included  here  and  the  interested  reader  is  encouraged  to  con¬ 
sult  reference  2  for  more  detail. 

14 

Consider  the  two  families  of  p-v  relations  depicted  in  Figure  7. 

Both  models  are  designed  so  as  to  allow  the  relative  amount  of  hardening 
and  softening  in  the  p-v  behavior  to  be  adjusted  by  varying  a  few  parame¬ 
ters.  From  the  extensive  study  conducted  [2],  we  draw  the  following  con¬ 
clusions  : 

(1)  relations  with  hardening  and  softening:  It  is  found  that  only 
a  limited  amount  of  hardening  is  permissible  for  a  solution  to  exist. 

To  demonstrate  this  finding,  consider  the  following  heuristic  argu¬ 
ment  in  connection  with  Figures  8a-e.  Assume  first  that  a  reasonable 
craze  profile  as  shown  in  Figure  8a  is  the  result  of  prescribing  a  p-v 
relation  depicted  in  Figure  8b.  Then,  without  considering  continuum 
mechanics,  one  deduces  from  8a  and  b  the  cohesive  force  distribution  in 
Figure  8c.  However,  on  mechanical  grounds  i.e.  using  equation  (2.19), 
this  p(x)  gives  rise  to  a  displacement  profile  v(x)  as  shown  in  Figure  8d 
which  is  not  equal  to  that  in  Figure  8a  and  exhibits  interpenetration  in 
the  craze  tip  region.  Moreover,  the  p-v  relation  deduced  then  from 

14.  Note  that  in  both  models,  the  fibril  stress  vanishes  for  displacements 
greater  than  a  critical  value  v^  ;  we  thus  deal  with  a  critical  crack  tip 
opening  displacement  (CTOD)  criterion  which  can  be  Justified  on  physical 
grounds  as  is  discussed  at  length  in  part  2  of  this  work  [1]. 

15.  Quatltative  results  are  given  in  reference  2. 
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Figures  8c  and  8d  and  shown  in  Figure  8e  bears  little  resemblance  to  the 
Initial  p-v  relation  of  Figure  8b.  Thus,  through  this  somewhat  roundabout 
argument,  we  conclude  that  when  the  prescribed  p-v  relation  hardens  by  an 
'appreciable'  amount  prior  to  softening  (as  in  Figure  8b)  and  our  solution 
algorithm  diverges  (i.e.  cannot  find  a  solution)  for  a  wide  range  of  ini¬ 
tial  guesses,  there  is  no  solution  for  the  p-v  relation  prescribed.  This 
line  of  argument  is  helpful  in  ascertaining  the  nonexistence  of  a  solution 
when  our  algorithm  falls  to  produce  one. 

(2)  Monotonical ly  softening  p-v  relations  (v  =  oi:  We  observe  that 

m 

1 6 

when  exceeds  certain  values,  only  the  trivial  solution  results,  oth¬ 
erwise  a  nontrivial  solution  exists.  A  heuristic  argument  similar  to  the 
one  presented  in  finding  (1)  above  is  given  in  reference  2. 

We  do  not  address  the  question  of  the  uniqueness  of  solutions 
rigorously.  As  Is  the  case  in  most  nonlinear  problems,  uniqueness  is  dif¬ 
ficult  to  establish.  Nevertheless,  based  on  the  extensive  experience 
acquired  during  the  course  of  this  study,  we  are  convinced  that  for  a 
given  craze  length  and  a^,  the  solution  found  by  using  Picard's  algorithm 
is  unique.  However  we  find  that  for  a  given  craze  length,  there  may  be 
two  different  values  of  <r^  and  two  corresponding  sets  of  p(x)  and  v(x) 
which  form  solutions:  The  solution  with  the  higher  corresponds  to  a 

full  craze  or  a  craze  with  a  very  short  internal  crack,  while  the  other 
one  with  the  lower  <s^  is  Invariably  the  solution  of  a  craze  with  a 
'longer'  central  crack.  These  aspects  are  discussed  In  more  detail  in 
references  1  and  2. 


16.  The  solution  corresponding  to  a  fully  closed  craze,  i.e.  v(x)  =  0  at  all 
points  or,  simply:  No  craze  exists. 


5.  CONCLUDING  REMARKS 


We  have  shown  that  the  problem  of  a  craze  (with  or  without  an  inter¬ 
nal  crack)  located  in  an  infinite  plane  under  craze-normal  loading  can  be 
solved  exactly  within  numerical  limitations.  The  formulation  and  solution 
scheme  are  simple  and  straightforward.  The  implemented  Picard's  algorithm 
is  effective  and  much  more  economical  than  Newton's  method.  The  formula¬ 
tion  can  also  be  used,  with  definite  advantage  over  the  presently  avail¬ 
able  schemes,  to  determine  stress  distribution  profiles  and  the  p-v  rela¬ 
tions  from  experimental  displacement  profiles.  Certain  characteristics  of 
the  fibril  force-displacement  (p-v)  relation  which  are  incompatible  with 
observed  craze  profiles  are  pointed  out.  The  methodology  developed  here 
is  applied  in  Part  2  of  this  study  [1]  in  simulating  craze  and  crack 
growth,  yielding  an  additional  set  of  informative  findings. 
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8.  FIGURE  CAPTIONS 


Figure  1.  Electron  micrographs  of  a  craze  (a)  and  its  tip  region  (b). 
(Courtesy  of  L.  Berger  and  E.J.  Kramer.) 


Figure  2.  (a)  Problem  geometry  definition 

(b)  Actual  craze  contour  with  primordial  craze  contour 

(c)  Linearized  geometry:  net  displacement  profile 

(d)  Model  for  craze  fibril  p-v  relation 


Figure  3.  Experimental  craze  contour  from  reference  7. 


Figure  4.  Stress  distribution  profile  derived  from  V(X)  in  Figure  3 
using  the  method  in  Section  3. 


Figure  5.  The  p-v  relation  derived  from  Figures  3  and  4. 

Figure  6.  Convergence  illustration  v  versus  x  (a®  =  0.05,  r  =  1.5). 

Iterations  #  2-8  fall  In  between  iterations  #  1  &  9. 


Figure  7.  Two  p-v  relations  used  in  studying  characteristics  of 
' admissible ' p-v  relations. 

(a)  piecewise-llnear  model 

(b)  cosine  model 


Figure  8.  Demonstration  of  the  nonexistence  of  solution  for  prescribed 
p-v  relations  such  as  the  one  in  figure  (b)  .  See  discussion 
in  Section  4.2. 


(b) 


Figure  1  Electron  microgrophs  of  a  craze  (a) 
and  its  tip  region  ( 0 ) . 

(Courtesy  of  L.  Berger  and  E.J.  Kramer.) 


(b) 


V(X)=  Net  displacement 
profile  w.r. t. 


the  linearized 
B.C.’s 


(d) 

(a)  Problem  geometry  definition 

(b)  Actual  craze  contour  with  primordial  craze  contour 

(c)  Linearized  geometry;  net  displacement  profile 

(d)  Model  for  craze  fibril  p- v  relation 


Figure  2. 


Figure  3.  Experimental  craze  contour  from  Lauterwasser  and  Kramer 
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ABSTRACT 


In  this  study  we  investigate  the  effects  of  nonlinear  fibril  behavior 
on  the  mechanics  of  craze  and  crack  growth.  The  effect  of  strain¬ 
softening  cohesive  material  on  crack  stability  is  of  particular  Interest 
and  is  examined  via  a  craze  and  crack  model  developed  in  reference  1  where 
the  formulation  and  solution  of  the  problem  are  discussed.^ 


In  this  second  part,  quasi-static  growth  of  a  craze  with  a  central 
crack  is  analyzed  for  different  nonlinear  force-displacement  (p-v)  rela¬ 
tions  for  the  craze  fibrils.  A  'critical  crack  tip  opening  displacement' 
(CTOD)  or  more  precisely,  'critical  fibril  extension'  is  employed  as  the 
criterion  for  fracture.  The  p-v  relation  is  further  assumed  to  be  invari¬ 
ant  with  respect  to  the  craze  and  crack  lengths.  The  results  are  compared 
with  the  Ougdale  model;  the  craze  zone  size  and  the  energy  dissipation 
rate  approach  asymptotic  values  in  the  limit  of  long  cracks. 


The  problem  of  craze  growth  from  a  precut  crack  under  increasing 
far-fleld  loading  is  then  studied.  In  the  case  where  the  p-v  relation  is 
monotonically  softening,  the  crack  can  start  to  grow  in  an  unstable  manner 
before  the  crack  tip  opening  displacement  reaches  its  critical  value. 


1.  Here  all  quantities  are  presented  in  their  nondlmensionalized  forms  as 
defined  in  reference  1.  The  unabridged  version  of  this  work  is  given  in 
reference  2. 


'j*  'd 


1.  INTRODUCTION 


In  reference  1  the  mathematical  and  experimental  treatment  of  craze 
formation  were  summarized.  We  developed  there  an  analytical  and  computa¬ 
tionally  efficient  scheme  to  obtain  the  stress  distribution  and  displace¬ 
ment  profiles  of  a  craze/crack  for  a  prescribed  craze  length  and  an  arbi¬ 
trarily  nonlinear  force-displacement  (p-v)  behavior  of  the  craze  fibrils. 
The  solution  scheme  developed  in  references  1  and  2  is  here  applied  to 
investigate  how  the  nonlinear  fibril  p-v  characteristics  affect  the 
mechanics  of  craze  and  crack  growth.  The  critical  crack  tip  opening  dis¬ 
placement  (CTOD)  is  employed  In  a  natural  way  as  a  fracture  criterion. 
The  applicability  of  this  CTOD  criterion  to  crazes  containing  cracks  is 
based  on  experimental  evidence  and  will  be  discussed  in  the  next  section. 
It  is  assumed  that  under  quasi-static  conditions  the  fibril  p-v  relation 
does  not  change  its  character  as  the  craze  and  crack  propagate.  This 
assumption  is  made  for  lack  of  experimental  data  and  because  of  an  inade¬ 
quate  understanding  of  the  fibril  behavior  at  the  present  time.  In  a 
later  section,  we  examine  the  problem  of  a  craze  initiating  from  an  exist¬ 
ing  (precut)  crack  with  special  interest  on  the  effect  of  'softening' 
fibril  behavior  on  crack  stability. 


2.  THE  TWO  BASIC  ASSUMPTIONS 

We  examine  first  the  'crack  tip  opening  displacement’  criterion 
(CTOD)  and  the  'invariance'  assumption  regarding  the  fibril  p-v  behavior 
as  the  craze  and  crack  grow  in  the  light  of  existing  experimental  results. 
Of  these  separate  issues  we  address  first 

2.1.  The  CTOD  Criterion. 

When  a  central  crack  is  present  within  a  craze,  the  crack  may  or  may 
not  propagate  as  the  far-fleld  load  is  Increased.  If  the  crack  does 
not  propagate,  it  is  clear  that  the  fibril  at  the  crack  tip  can  sustain 
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additional  growth  (lengthening)  though  not  necessarily  higher  stress  or 
strain.  This  is  the  consequence  of  two  physical  mechanisms  that  operate 
simultaneously,  namely,  surface  drawing  and  fibril  creep.  If  surface 
drawing  alone  operates,  a  fibril  will  never  break  and  the  crack  will  never 
advance.  This  possibility  is  contrary  to  experimental  observations  since 
cracks  do  grow  within  crazes.  As  can  be  explained  in  terms  of  craze  and 
fibril  microstructurcs ,  a  crack  that  propagates  quasi-statically  through  a 
craze  either  runs  right  along  the  midrib  where  the  fibrils  are  weakest  or 
propagates  in  a  'patch'  or  'mackerel'  pattern  [3,4].  In  the  first  situa¬ 
tion,  the  heightened  stress  at  the  crack  tip  causes  the  previously 
unloaded  midrib  to  stretch  until  the  local  extension  ratio  exceeds  the 
maximum  sustainable  by  the  bundles  of  polymeric  molecule  chains  that  con¬ 
stitute  the  fibril,  thus  leading  to  fracture.  In  the  latter  event,  usu¬ 
ally  observed  at  higher  crack  velocities  [4],  the  high  crack  tip  stress 
causes  fibrils  at  the  craze-polymer  interface  to  be  drawn  at  high  local 
strains  [3].  One  observes  then  that  the  crack,  in  this  case,  prefers  to 
run  along  either  Interface  where  the  fibrils  are  locally  weaker  than  at 
the  midrib,  creating  the  familiar  mackerel  pattern  [3,4].  In  either  case, 
the  above  situations  support  the  use  of  a  critical  crack  tip  opening  dis¬ 
placement  (or  a  critical  fibril  extension  criterion).  This  is  also  con¬ 
firmed  by  experiments  performed  by  Doll  et  al  [5,6]  who  found  that  the 
CTOD  is  relatively  constant  for  a  wide  range  of  crack  velocities.  Specif¬ 
ically,  these  authors  showed  that  the  maximum  craze  widths  at  the  crack 
tip  (critical  CTOD  or  v  in  our  notation),  measured  in  PMMA  compact  ten- 
sion  specimens  for  crack  speeds  ranging  from  10  to  20  mm/sec,  are  essen- 
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tially  constant  with  a  standard  deviation  of  about  8S5.  Their  findings 

2.  Surface  drawing  is  the  process  by  which  bulk  polymer  is  converted  (drawn) 
into  fibrillar  craze  material.  Fibril  creep  is  simply  the  stretching  of 
an  existing  fibril  without  drawing  additional  polymer  material  into  the 
craze  from  the  bulk. 

3.  In  private  communication.  Professor  Kramer  indicated  that  under  a  low 
crack  speed  of  about  1  micron/sec  the  critical  CTOD  increases  substan¬ 
tially.  This  casts  some  doubt  on  the  validity  of  the  number  10~®  mm/sec 
given  in  references  5  and  6.  However,  for  the  speed  range  between  10~^ 
and  20  mm/sec,  the  constancy  of  the  CTOD  (for  PMMA)  appears  acceptable. 
See  the  introduction  of  reference  13  for  further  discussion  on  this  point. 


suggest  that  the  fibrils  spanning  the  tip  of  the  crack  can  ojily  sustain  a 
limited,  relative,  displacement  of  their  ends.  An  interpretation  of  the 
experimental  result  in  terms  of  molecular  behavior  is  given  in  reference 
6;  such  discussions  are  beyond  the  scope  of  the  present  work.  For  our 
purpose,  it  suffices  to  assume  that  the  surface  drawing  mechanism  is 
somehow  constrained  at  the  crack  tip  and  creep  invariably  gives  rise  to 
fibril  breakage  at  a  roughly  constant  CTOD.'*  We  assume  in  this  study, 
therefore,  that  for  a  given  polymer,  a  critical  CTOD  exists  and  is  con¬ 
stant  . 

2.2.  Comments  on  the  Invariance  of  the  p-v  Relation. 

To  study  craze  and/or  crack  growth,  one  needs  to  know  whether  the  p-v 
relation  changes  as  propagation  proceeds  and  if  so,  how.  The  question  can 
only  be  answered  by  experiments.  Since  no  experimental  data  on  this  par¬ 
ticular  topic  is  available,  we  feel  free  to  assume  further  that  the  p-v 
relation  is  'invariant'  with  respect  to  both  the  craze  and  the  crack 
lengths  for  quasi-static  propagation  under  a  'controlled'  environment.  As 
it  turns  out,  on  the  basis  of  the  subsequent  studies,  this  assumption  can 
only  be  valid  over  a  limited  range  of  craze  and  crack  lengths  (see  Section 
3). 

3.  NUMERICAL  SIMULATIONS  OF  CRAZE  AND  CRACK  PROPAGATION 

We  proceed  now  to  apply  the  solution  scheme  of  reference  1  to  study 
how  the  nonlinear  fibril  (cohesive)  p-v  characteristics  affect  the  mechan- 
Ics  of  craze  and  crack  growth.  The  geometry  of  the  problem  is  depicted 

4.  It  is  worthwhile  to  note  that  cracks  are  observed  to  form  in  'full'  crazes 
at  the  polymer/craze  interface  and  not  at  the  midrib.  The  locations  of 
the  first  formed  crack  in  a  full  craze  are  also  not  necessarily  at  the 
center  part  of  the  craze  where  the  fibrils  are  longest. 

5.  Temperature,  pressure,  moisture  content,  chemical  environment  (such  as  the 
presence  of  alcohol  or  any  chemical  agents  possibly  affecting  fibril  draw¬ 
ing  and  creep)  are  all  kept  unchanged. 


in  Figure  1:  The  craze  length  is  denoted  by  c,  the  crack  length  by  'a'  and 

the  vertical  boundary  displacement  of  the  craze  by  v.  The  remotely 
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applied  stress  is  and  the  fibril  cohesive  force  is  represented  by  p. 

Given  a  craze  length  and  a  p-v  relation  for  the  fibrils,  the  problem  is  to 
find  p(x)  ,v(x)  and  <s^  such  that  both  equilibrium  and  the  smooth  closure 
(Barenblatt ' s)  condition  are  satisfied  simultaneously.®  The  solution 
scheme  is  presented  in  references  1  and  2. 

Quasi -static  craze  and  crack  growth  is  simulated  by  varying  the  craze 
length  c  continuously.  Based  on  the  discussion  in  Section  2,  we  assume 
the  fibril  force-displacement  relation  of  the  form 

p(x)  =  p(v(x)]  (3.1) 

The  effect  of  several  p-v  relations  will  be  examined.  We  start  by 
first  reviewing  the  Barenblatt-Dugdale  model  [7,8]  shown  in  Figure  2. 

3.1.  The  Barenblatt-Dugdale  Model. 

This  simplest  model  assumes  constant  cohesive  (fibril)  stress 

throughout  the  cohesive  (craze)  zone  as  denoted  by  p^  in  Figure  2.  The 

problem  has  been  analyzed  by  Rice  [9]  and  Goodier  and  Field  [10];  we  sum- 

g 

marize  below  the  pertinent  results  for  reference.  They  are  obtained 
directly  from  the  equations  presented  In  Section  2  of  reference  1.  The 

6.  We  consider  here  the  planar  problem  i.e.  plane  stress  or  plane  strain. 

7.  The  mathematical  modelling  of  the  problem  has  been  described  in  references 
1  and  2. 

8.  In  this  paper,  the  (net)  stress  Intensity  factor  is  always  zero  and  thus 

does  not  enter  the  discussion  as  a  pertinent  parameter.  The  external 
loading  is  completely  characterized  by  which  is  related  to  the  'load¬ 
ing'  stress  intensity  factor  by  \|nc  .  Note  that  the  stress 

states  at  both  the  crack  tip  and  the  craze  tip  are  determined  by  the  non¬ 
linearity  of  the  fibril  p-v  relation  rather  than  the  'usual'  K-field. 

9  All  equations  and  parameters  are  in  dimensionless  forms. 
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Barenblatt  condition  requires 


2  -l/a> 

=  z  P„cos  (-) 
o"  n  m  c 


(3.2) 


and  the  displacement  v(x)  is 


r 


v(x)  =  !  alog  i  ^ 


2  2  , 
c  -X  -t-X 


2  2  1„2  „2 


I 


c^-a^  ,  .  ,  x\  lc‘--a‘--a\lc‘--x^  ,  1 

'  I  +  xlog  I  - - -  I 


I 


.  2  2,22 

\  c  -X  -\  c  -a 


,22,22  (3.3) 

x\  c  -a  +a\  c  -X  J  '  ' 


At  x=a,  V  equals  the  critical  CTOD,  thus  applying  L'Hospital's  rule, 


one  obtains 


2P„ 


v(a)  =  ^  a  iog(^) 

c  n  a 


(3.4) 


Let  «  denote  the  size  of  the  craze  (yield)  zone,  i.e. 


«  H  c  -  a 


(3.5) 


Further,  let  7  be  the  'fibril  fracture  energy'  (or  the  'surface  energy') 
which  Is  the  work  expended  In  pulling  a  fibril  from  v  =  0  to  v  =  v^.  For 
the  Barenblatt-Dugdale  model,  we  have  simply 


y  =  P_v, 


(3.6) 


By  analyzing  the  data  for  six  polymers  in  references  3  and  11,  we 


find  that  the  dimensionless  parameters  p  and  v„  have  values  in  the  ranges 

®  c 


of 


0.01  <  p  ^0.1 

m 


0.001  <  V  <;  0.01 

c 


We  shall  use  these  numbers  later  as  guides  for  bounding  the  various  param¬ 
eters  in  the  p-v  relations  to  be  studied  so  that  the  results  of  our 


.•■-j 


analyses  will  be  quantitatively  meaningful. 


We  examine  next  the  energy  release  rate  of  the  system.  Let  U 
represent  the  elastic  energy  of  the  system.  The  energy  release  rate  G  is 
defined  as  in  terms  of  the  potential  energy  n  as 


G  = 


in 

da 


(3.7) 


Consider  the  fixed  grip  case  so  that  no  external  work  is  done  during  the 
course  of  crack  propagation.  Then  (3.7)  reduces  to 


G  = 


ay 

da 


(3.8) 


Now  U  Is  given  by  the  elastic  energy  stored  initially  in  the  plate 
(planar  body),  U^,  minus  the  energy  which  is  required  to  form  a  crack 
of  length  2a  and  two  craze  zones  of  length  m  each.  We  have  then 


U  =*  Up  -  Up  (3.9) 

=  i  (3.10) 

P  4  oo 


where  we  let  d  denote  the  dimension  of  this  'infinite'  square  plate  (i.e. 
d  »  c).^®  U  is  given  by 

V 


Uc  =  4  [  -ya  ^  ;  p^v(x)dx  ]  ^3 

where  the  factor  of  4  accounts  for  the  four  quadrants.  Substituting  v(x) 
from  (3.3)  and  making  use  of  (3.2)  and  (3.4),  one  obtains  after  some  alge¬ 
bra 


10.  Note  that  the  factor  1/4  in  place  of  the  normally  expected  factor  of  1/2 
is  a  result  of  the  particular  non-dimensional ization  scheme  employed. 


-  9  - 


=  4a[  7  -t-  y  {  \ic^-a^  -  2v^  } 


(3,12) 


Using  (3.9)  and  (3.10)  in  (3.8)  yields 


,  „  a<T  3U 

P  _  _1  -  ,2  .  _ ^  _ !2 

2  “  Ba  da 


(3.13) 


At  this  point  an  important  question  enters:  What  is  the  criterion 
for  crack  propagation  (fracture)  ? 


Goodier  and  Field  [10]  assumed  that  the  crack  grows  in  a  self-similar 
manner,  such  that 


ac  _  c 
da  a 


(3.14) 


This  condition  is  the  consequence  of  assuming  a  to  be  constant  as  the 

oo 

crack  and  craze  grows.  Under  these  conditions,  -r —  vanishes  and  G 

da 

^^c 

becomes,  in  view  of  (3.13),  which  is  the  'plastic'  (cohesive)  work 

rate.  However,  the  crack  tip  opening  displacement  v(a)  =  v^  as  given  by 
(3.4)  becomes  a  linear  function  of  a.  Thus  the  CTOD  in  the  model 
presented  by  Goodier  and  Field  model  increases  as  the  crack  grows.  This 
behavior  is  inconsistent  with  experimental  observations  as  discussed  in 
Section  2. 


In  the  present  model,  the  critical  CTOD,  v^,  is  kept  constant;  this 
condition  in  turn  requires  that  must  decrease  as  'a'  increases  (see 
Figure  7)  in  order  to  maintain  the  quasi-static  condition,  i.e. 


(3.15) 


Viewing  (3.13)  in  this  light,  it  is  evident  that  G  becomes  unbounded  if  d 

da  3U 

oo  I*' 

is  unbounded,  no  matter  how  small  a^{-r — )  may  be,  since  is  always  fin- 
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Ite. 


Therefore,  a  crack,  the  tip  displacement  of  which  has  reached  the 
critical  value,  is  always  unstable  in  this  geometry. To  this  result  we 
only  need  to  add  the  reminder  that  for  a  clamped  infinite  (very  large) 
plate,  a  small  enlargement  of  an  Internal  crack  would  hardly  reduce  the 
far-field  <T^  sufficiently  to  keep  the  crack  and  craze  growing  in  a  quasi¬ 
static  manner.  Therefore  once  a  crack  (with  or  without  a  cohesive  zone) 
starts  growing,  and  the  fracture  criterion  such  as  a  critical  CTOD  is 
satisfied,  catastrophic  failure  ensues.  We  shall  see  later  that  the  con¬ 
verse  is  not  true. 


The  plastic  work  rate,  >  calculated  using  the  CTOD  criterion 
differs  markedly  from  the  result  obtained  by  Goodier  and  Field  under  the 
assumption  of  self-similar  growth  (see  [2]  and  [10]  for  details).  Based 
on  the  present  criterion,  the  plastic  work  rate  increases  as  'a'  increases 
and  approaches  the  constant  value  47  rapidly.  Physically,  for  a  short 
crack,  one  has  a  relatively  large  craze  zone  w  ;  to  advance  the  crack  by 
an  amount  Aa,  very  little  energy  is  needed  since  c  increases  more  slowly 
than  'a'  (l.e.  the  craze  zone  w  shrinks).  For  longer  cracks,  the  size  of 
the  craze  zone  tends  toward  a  constant  value  and  thus  the  craze  zone  sim¬ 
ply  translates  as  the  crack  propagates.  In  this  latter  case,  the  dissi¬ 
pated  energy  is  thus  47  per  unit  crack  advance.  This  is  demonstrated 
later  in  Figure  8  together  with  results  for  other  nonlinear  p-v  models. 

3.2.  General  p-v  Models. 

Next  we  examine  six  different  p-v  relations  labeled  #1  through  #6  in 
Figure  3  and  compare  the  results  to  the  Dugdale  model  labeled  Model 

#6  approximates  the  p-v  behavior  of  craze  fibrils  as  it  is  known  today 

11.  If  the  tip  displacement  is  less  than  v^ ,  the  crack  does  not  grow  as  the 
tip  fibril  has  not  yet  been  broken.  For  softening  p-v  relations,  however, 
instability  may  occur  before  v^  is  reached  (see  Section  4). 
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[11,12].  As  a  check,  we  employ  this  model  to  simulate  craze  and  crack 
growth.  The  resulting  craze  contours  and  stress  distribution  profiles  are 
depicted  in  Figures  4  and  5  respectively.  They  do  closely  resemble  the 
experimental  v(x)  and  p(x)  obtained  by  Kramer  and  his  coworkers  [11,12]. 

To  compare  the  results  for  all  these  different  p-v  relations,  it  is 

desirable  to  use  the  Dugdale  model  as  a  benchmark  (here  we  use  the  model 

labeled  'D'  in  Figure  3  with  p  =  0.025  and  v  =  0.008  ).  All  the  p-v 

m  c 

models  are  designed  so  that  the  area  under  each  p-v  curve,  y  ,  is  constant 
and  equal  2x10  the  same  as  that  for  the  Dugdale  model  employed. 


The  numerical  procedures  have  been  presented  in  Section  2.3  of  refer¬ 
ence  1.  The  craze  length  is  varied  from  c  =  0.0  to  c  =  10.0,  the  incre¬ 
ment  Ac  ranges  from  0.01  to  0.2.  Smaller  increments  for  Ac  are  necessary 
to  capture  the  behavior  of  the  system  for  short  crack  lengths,  particu- 

^^C  12 

larly  in  calculating  the  rate  of  plastic  work,  During  the  course  of 

these  craze  growth  simulations,  it  was  found  that  convergence  is  generally 


attained  in  fewer  iterations  if  Vj^(c)  is  used  as  the  Initial  guess  of 
v^(c+Ac).  (This  is  not  the  case  for  p-v  model  #6,  however.)  In  order  not 
to  confuse  the  presentation  of  the  results,  the  trivial  solution  which 
occurs  at  short  craze  lengths  and  corresponds  to  a  closed  craze  (v(x)  =  0) 
is  not  Included.  Also,  for  the  Barenblatt-Dugdale  model  'D',  results  that 
are  'nonphysical'  in  the  sense  of  a  crack  growing  in  a  shrinking  craze, 
are  omitted.  Interested  readers  are  refered  to  reference  2  for  more 
details . 


In  Figure  6,  the  craze  zone  u  Is  plotted  as  a  function  of  the  crack 
length  'a*.  The  craze  zone  is  seen  to  decrease  monotonically  and  to  tend 

1  3 

toward  constant  asymptotic  values  as  'a'  Increases  for  all  p-v  models. 


12.  is  obtained  by  integrating  numerically  ^ 

U  =  4[  ya  *■  f  ev(v(x))dx]  where  ev(v)  =  /  p(v')dv' 

^  a  o 

13.  For  the  Barenblatt-Dugdale  model,  this  can  be  analytically  shown  to  be  the 


Thus  in  the  limit  of  long  cracks,  the  craze  zone  is  simply  shifted  as  the 
crack  and  craze  grow.  Note  that  w  is  primarily  dependent  on  but  also 
depends  weakly  on  the  character  of  the  p-v  relation. 

In  Figure  7,  the  far-fleld  stress  <s  normalized  by  p  =  0.025  is 
plotted  versus  c  for  the  various  p-v  models.  The  remotely  applied  stress 
decreases  monotonlcally  as  the  craze/crack  grows  which  is  necessary  to 
maintain  quasi-static  growth.  This  confirms  the  conclusion  reached  in  the 
previous  section  from  the  energy  release  rate  analysis  that  once  the  crit¬ 
ical  CTOD  is  reached  the  crack  will  grow  unstably  if  the  remotely  applied 
stress  is  maintained  constant.  An  Interesting  feature  is  exhibited  by  the 
dash-dot  curve  which  represents  a  'full  craze'  (a  =  0)  solution  branch  and 
is  present  only  in  the  p-v  model  #6  (among  all  models  studied  here).  This 
unusual  behavior  merits  further  elaboration: 

One  finds  that  for  a  craze  length  c  <  0.8,  the  algorithm  converges  to 
the  trivial  solution  in  which  the  craze  is  closed  (v(x)=0)  at  all  points 
(see  Figure  13).  For  c  greater  than  0.8,  there  exists  a  'full  craze' 
solution  regime  up  to  c  =  2.2,  At  c  =  2.2,  the  solution  shows  that  v(o) 
is  approximately  0.006  and  thus  p(x)  has  a  minimum  at  v(x)  =  0.005  (see 
Figure  3).  For  c  larger  than  2.2,  there  is  no  full  craze  solution  (the 
algorithm  'jumps'  and  converges  to  solutions  with  cracks  i.e.  a  5*  0). 
This  behavior  can  again  be  explained  using  a  heuristic  argument,  see  Sec¬ 
tion  3.2  of  reference  2.  Note  that  for  1.87  <  c  <  2.2  ,  the  solutions  are 
not  unique,  i.e.  for  a  given  c  in  this  range,  we  either  have  a  full  craze 
solution  at  a  higher  a  or  a  crack/craze  solution  at  a  lower  a  . 

OO  00 

The  above  result,  which  indicates  the  absence  of  solutions  for  the 
prescribed  p-v  relation  (#6)  over  some  ranges  of  craze  and/or  crack 
lengths  may  seem,  at  first  sight,  contradictory  to  our  physical  intuition: 
One  would  think  that  crazes  and/or  cracks  of  different  sizes  should  exist 
in  a  continuous  size  'spectrum'.  This  seeming  contradiction  is  resolved 
by  recognizing  that,  most  probably,  the  p-v  relation  is  not  'invariant' 
with  respect  to  all  values  of  the  crack  and  craze  lengths  as  initially 


13 


assumed. More  elaborate  modelling  of  the  fibril  p-v  characteristic  is 
therefore  needed;  see  Section  5  for  further  discussion. 


In  Figure  8,  the  rate  of  'plastic'  work. 


da 


is  plotted  as  a  func¬ 


tion  of  the  crack  length  'a'.  It  approaches  the  asymptotic  limit  47 


fairly  quickly  for  all  p-v  models  studied. 
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4.  CRAZE  GROWTH  INITIATING  FROM  A  PRECUT  CRACK 

In  this  section,  we  consider  the  problem  of  a  craze  growing  from  a 
precut  slit  as  shown  in  Figure  9.  Under  the  assumption  of  a  CTOD  cri¬ 
terion,  one  would  expect  the  crack  to  start  advancing  when  the  remotely 

applied  stress  <t  is  high  enough  to  make  the  crack  tip  displacement  v(a) 
00 

reach  the  value  v  .  Thereafter  the  crack  and  craze  would  grow  unstably  as 
c 

discussed  In  the  previous  section. 

However,  it  is  conceivable  that  for  certain  p-v  relations,  a  may 

00  * 

reach  its  maximum  value  before  the  crack  tip  opening  displacement  reaches 
its  critical  value  v^  ,  thus  resulting  in  instability  before  the  CTOD  cri¬ 
terion  Is  satisfied.  We  shall  show  that  this  is  Indeed  the  case  for  p-v 
relations  that  soften  monotonlcally  up  to  fracture  (fibril  breakage).  To 
start,  let  us  first  study  the  Barenblatt-Dugdale  model  for  later  refer¬ 
ence. 


14.  Note  that  this  conclusion  is  based  solely  on  mechanics  analysis. 

15.  The  curve  representing  p-v  model  #6  has  been  smoothed  out  for  a  >  2.5. 
For  this  particular  p-v  relation,  to  get  accurate  v^  in  the  craze  zone 
(whose  size  is  =  0.8  for  a  >  2.5),  finer  discretization  is  required  to 
capture  the  steep  rise  and  fall  of  the  stress  distribution  in  the  craze 
zone . 
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4.1.  The  Barenblatt-Dugdale  Model . 

We  start  with  a  precut  crack  of  length  2a^  Initially  the  craze 

length  c  =  a  ,  l.e.,  w  =  0.  We  then  increase  a  gradually  until  v(a) 

o 

reaches  v^.  By  equation  (3.4)  when  v(a)  =  v^,  the  'critical'  craze  length 

is 

cr 

TIV 

‘^cr  ^  2p  a  (4.1) 

m  o 

Hence  for  a  s  c  <  c  ,  the  remotely  applied  stress  ff  can  be  obtained 
o  cr 

from  (3.2),  i.e. 


<r 

00 


-1 

cos 


a 

o 


5  c  <  c 


cr 


(4.2) 


Based  on  (4.1)  and  (4.2),  we  can  plot  the  applied  stress  as  a  function 
of  c  for  constant  a  =■  a^  as  shown  in  Figure  10.  Four  initial  crack 
lengths  are  used  (a^=  0.6,  0.9,  1.2,  and  1.5).  The  curve  representing 
quasi-static  crack  and  craze  growth  referred  to  subsequently  as  the  'mas¬ 
ter*  curve)  and  on  which  the  CTOD  always  attains  its  critical  value  is 
also  identified.  On  the  'master'  curve,  the  craze  and  crack  growth  is 
unstable  in  the  sense  that  must  be  reduced  continually  in  order  to 
maintain  quasi-static  propagation. 


Along  the  curves  on  which  a  =  a^,  the  CTOD  v(a)  gradually  grows  as 
the  remotely  applied  stress  is  increased  e.g.  from  point  A  to  B;  it 

then  reaches  its  critical  value  v^  at  B.  Thus  for  the  Barenblatt-Dugdale 
model,  <r  Increases  monotonicai ly  in  the  process  and  instability  can  occur 

oo 

only  when  the  'master'  curve  is  reached. 


4.2.  General  p-v  Models. 


Next  let  us  consider  the  same  problem  as  in  Section  (4.1)  using  p-v 
relations  #1  through  #6  of  Figure  3.  In  contrast  to  the  Barenblatt- 
Dugdale  model,  we  find  that  for  p-v  relations  that  soften  monotonical ly  up 
to  the  point  of  fracture  (models  #1  through  #5)  instability  may  occur 
before  the  critical  CTOD  is  reached  at  the  crack  tip.  Briefly,  the  numeri¬ 
cal  analysis  proceeds  in  the  following  sequence: 

1.  Start  with  a  value  a^  which  corresponds  to  a  c^^.  that  lies  on  the 
master  curve. 

li.  Increase  c  gradually  from  a^  to  c^^  in  such  way  that  the  crack  tip 
falls  exactly  on  a  discretized  point  x^  .  This  avoids  unnecessary 
Interpolation . 

iii.  Find  the  solution,  pre-specifying  that  p(x)  vanish  for  x  5  a  The 

o 

solution  would  automatically  have  v(a)  <  v  as  long  as  c  <  c 

c  cr 

For  Illustrative  purposes  we  present  first  the  result  for  model  #5  as 

being  representative  of  the  class  of  monotonically  'softening'  p-v  models 

(#1  through  95).  More  complete  results  are  given  in  [2].  In  Figures  11 

and  12  (enlargement  of  a  part  of  Figure  11),  one  secs  that,  for  a  =03 

o 

to  1.5,  instabilities  occur  before  the  'master'  curve  is  reached. 

Consider  a  stress-controlled  experiment.  As  the  far-fleld  stress  a 

oo 

is  Increased  from  zero  to  point  B  (Figure  12),  the  CTOD,  l.e.  v(a), 

approaches  the  critical  CTOD,  v  .  At  B,  any  further  increase  in  a  will 

C  ox 

cause  the  crack  and  craze  to  grow  unstably  since  quasi-static  growth  will 
no  longer  be  maintained  at  that  stress  level.  If  one  reduces  a  instead 

OO  * 

keeping  in  mind  that  the  craze  length  c  cannot  decrease  while  a  =  a 

0 

remains  constant,  point  C  could  then  be  reached  quasi-statlcally .  There¬ 
fore,  at  point  B,  any  small  positive  perturbation  in  <J  or  in  a  would 

OO  O 

cause  the  crack  and  craze  to  grow  in  an  unstable  manner.  Note  that  the 
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solution  cannot  'jump'  spontaneously  from  point  A  to  point  C  because  more 
work  has  to  be  done  to  further  develop  to  the  craze  zone  and  this  is 
achieved  by  following  the  path  ABC  as  discussed. 

To  explore  this  instability  further,  consider  perturbing  the  system 

at  point  A  by  cutting  fibrils  at  the  crack  tip  making  a  increase  inflni- 

o 

tesimally.  Does  this  cause  catastrophic  failure?  As  depicted  in  Figure 
12,  the  solution  will  then  move  from  point  A  to  a  nearby  point  D  which 
corresponds  to  a^  =  0.3^  and  stop.  Thus  the  craze  and  crack  at  point  A  is 
stable  since  a  small  perturbation  in  a^  does  not  lead  to  total  failure. 
Note  in  passing  that  the  unstable  behavior  is  more  pronounced  for  the 
shorter  initial  crack  lengths. 

We  turn  next  to  the  result  for  p-v  model  #6  which  rehardens  after  the 
initial  softening  and  fails  abruptly  as  in  the  Barenblatt-Dugdale  model 
(see  Figure  13).  One  sees  that  the  unstable  character  observed  in  soften¬ 
ing  p-v  models  is  clearly  absent  which  closely  resembles  the  result  for 
the  Barenblatt-Dugdale  model  shown  in  Figure  10.  These  results  indicate 
that  in  cases  where  craze  fibrils  fail  in  an  abrupt  manner  (as  opposed  to 
gradual  softening  to  rupture),  instability  does  not  take  place  until  the 

CTOD  reaches  its  critical  value  v  . 

c 

5.  CONCLUDING  REMARKS 

In  this  second  part  of  the  study,  we  have  examined  the  quasi-static 
growth  of  cracks  and  crazes  by  considering  a  number  of  nonlinear  fibril 
p-v  models  along  with  the  Barenblatt-Dugdale  model.  The  energy  dissipa¬ 
tion  rate  consistent  with  the  critical  CTOD  criterion  is  established  for 
those  materials.  We  also  considered  the  case  of  craze  growth  starting 
from  a  precut  crack;  under  these  circumstances  instabilities  are  observed 
to  occur  before  the  critical  CTOD  is  reached  at  the  crack  tip  for  p-v 
models  that  soften  monotonlcal ly  to  fracture. 


As  discussed  earlier,  not  enough  is  known  about  the  p-v  characteris¬ 
tics  of  the  craze  fibrils.  Further  investigation  into  the  interplay 
between  the  surface  drawing  and  creep  mechanics  should  shed  some  light  on 
how  the  p-v  behavior  of  fibrils  changes  as  the  craze  and  crack  grow;  a 
first  attempt  along  this  line  has  been  undertaken  by  Kramer  and  Hart  [13]. 
Experiments  such  as  those  performed  by  Kramer  et  al.  [3,11,12]  need  to  be 
carried  out  for  a  wider  range  of  craze  and  crack  lengths  to  determine  to 
what  extent  the  fibril  p-v  relation  is  invariant  with  respect  to  c  and 
'a'  . 
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1.  Problem  definitions 


2.  The  Barenblatt-Dugdale  model. 


3.  p-v  models  used  in  the  nonlinear  craze  and  crack  propagation 
studies  in  Sections  3  and  4.  All  models  have  the  same  7  . 


4.  Craze  and  crack  propagation  simulations  for  p-v  model  #6:  Dis¬ 
placement  profiles. 


5.  Craze  and  crack  propagation  simulations  for  p-v  model  #6: 
Stress  distribution  profiles. 

6.  The  craze  zone  u  as  a  function  of  'a'  for  various  nonlinear 
p-v  models  (p-v  models  as  indicated). 


7.  —  as  a  function  of  c  for  various  nonlinear  p-v  models  (p-v 
models  as  indicated). 

8.  The  rate  of  plastic  work  dissipation  as  a  function  of  'a' 
for  various  nonlinear  p-v  models  (p-v  models  as  Indicated). 


9.  Craze  growth  initiating  from  a  precut  crack 

(a)  Initial  precut  crack,  no  far-field  loading. 

(b)  Increasing  far-field  loading,  no  crack  growth. 

(c)  Critical  CTOD  reached,  crack  growth  imminent. 

10.  as  a  function  of  c  for  a  =  0.6,  0.9,  1.2,  and  1.5  for  the 
Barenblatt-Dugdale  model  'D'. 

11.  —  as  a  function  of  c  for  =  0.3,  0.6,  0.9,  1.2,  and  1.5  for 
p-v  model  #5. 


Figure  12.  —  as  a  function  of  c  for  a^  =  0.3,  0.6,  and  0.9  for  p-v  model 

#§  (enlargement  of  a  part  of  Figure  11). 

Figure  13.  —  as  a  function  of  c  for  a  =  i.o,  1.2,  1.4,  1.6,  and  1.8  for 
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p-v  model  #6. 


tfatlon  simulations  fcir  p-v  model  #6;  Dis- 


The  rate  of  plastic  work  dissipation  as  a  function  of  'a 

for  various  nonlinear  p-v  models  (p“V  models  as  indicated)  . 
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Figure  11  —  as  a  fum  tion  of  i:  for 
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